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Evolution of the ESR absorption in a strong-leg spin ladder magnet (C 7 HioN 2 ) 2 CuBr 4 (abbrevi¬ 
ated as DIMPY) is studied from 300 K to 400 mK. Temperature dependence of the ESR relaxation 
follows a staircase of crossovers between different relaxation regimes. We ague that the main mecha¬ 
nism of ESR line broadening in DIMPY is uniform Dzyaloshinskii-Moria interaction (|D| = 0.20 K) 
with an effective longitudinal component along an exchange bond of Cu ions within the legs resulting 
from the low crystal symmetry of DIMPY and nontrivial orbital ordering. The same Dzyaloshinskii- 
Moriya interaction results in the lifting of the triplet excitation degeneracy, revealed through the 
weak splitting of the ESR absorption at low temperatures. 

PACS numbers: 75.10.Kt, 76.30.-v 
Keywords: low-dimensional magnet 


I. INTRODUCTION. 

Low-dimensional magnets are actively studied during 
last decades both theoretically and experimentally. Spin- 
ladder is one of the simplest models of the field, that is 
just one step more complicated than the Heisenberg spin 
chain, the keystone of the low-dimensional magnetism. 
Such a system consists of two chains forming the ’’legs” 
of the spin ladder, which are coupled by simple interchain 
coupling forming ’’rungs” of the ladder. The Hamiltonian 
of the single spin ladder with the equivalent positions 
along the ladder is 

T~L = Jleg + S2,iS2y+l) + 

i 

Jrung ^ ^ -|- '}~i anis (i) 

i 

it includes Heisenberg exchange couplings J/ eg and J rU ng, 
Zeeman interaction (with usually anisotropic g-tensor) 
and weak anisotropic spin-spin interactions TL a nis- 
Independently on the ratio between the J\ eg and J rung , 
the excitation spectrum of the spin ladder is gapped, 
ground state is non-magnetic and excited states are 
S = 1 quasiparticles A However, most of the experimen¬ 
tally available examples of the spin ladder systems are 
the so-called strong-rung ladders with the dominating 
in-rung interaction J rimg . Strong-leg ladders remains a 
rarity in this family. Additional complication of the real 
systems is a presence of the anisotropic spin-spin inter¬ 
actions breaking the ideal symmetry of the Heisenberg 
model. Such interactions limit excitations lifetime (to 
the point of total damping in some extreme cases^) and 
could lift degeneracy of the 5=1 states. Thus, esti¬ 
mation of such interactions strength and, ideally, search 
for the systems with negligible anisotropic part of the 


Hamiltonian is an important quest when comparing real 
magnets with model predictions. 

Adequate accounting for the effect of anisotropic in¬ 
teractions in a spin-gap magnet is also a challenge. This 
problem was addressed in a ID field theory models via 
bosonio^ and fermionio^ approaches and in an indepen¬ 
dently developed macroscopic models However reliable 
microscopic models remains a rarity (see, e.g. Ref. [b): 
most of the real spin-gap magnets have a complicated 
network of the exchange couplings allowing far too nu¬ 
merous possibilities of the anisotropic interactions pa¬ 
rameters. The adequate microscopic approaches are of 
particular interest in connection with a particular case of 
the effect of a uniform Dzyaloshinskii-Moria interaction 
on the properties of a quantum magnet^ 

Recently found organometallic compound 

(C 7 Hi 0 N 2 ) 2 CuBr 4 , abbreviated DIMPY for short, 
is an example of the strong-leg ladder with very weak 
anisotropic interactions —li Presence of the energy gap 
in the excitation spectrum was revealed by magnetic 
susceptibility^, specific heal£ and magnetizatio n 11 : 12 
bulk measurements as well as by inelastic neutron 
scatteririgj2il£ The energy gap was founc£ to be 
0.33 meV, it can be closed by the magnetic field 
OoHci = 3.0 T, while the saturation field is much 
higher— pL Q H sat « 30 T. The values of the exchange 
constants were determined from the DMRG fit of the 
measured inelastic neutron scattering spectra— and were 
found to be Ji eg = 1.42 meV and J rU ng = 0.82 meV. 
The magnetic field induced ordering is observed at very 
low temperatures ( T^ nax ' > sa 300 mK at fj, 0 H ~ 15 T) — 
Electron spin resonance (ESR) spectroscopy is a pow¬ 
erful tool to probe for the weak anisotropic interactions 
in the magnetic systems. Inelastic neutron scattering 
experiments^ have shown that DIMPY is an almost per- 
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feet realisation of the Heisenberg spin ladder. ESR tech¬ 
nique allows much higher energy resolution (routinely 
resolved ESR linewidth of 100 Oe corresponds approxi¬ 
mately to the energy resolution of 1 /reV) and thus allows 
to probe possible effects of anisotropic interactions with 
high accuracy. 

In the present manuscript we report results of the ESR 
study of low-energy spin dynamics in DIMPY in the tem¬ 
perature range from 400 mK to 300 K. We observe angu¬ 
lar and temperature dependences of the ESR linewidth 
at high temperatures which can be described as an effect 
of the uniform Dzyaloshinskii-Moria (DM) interaction, 
which is allowed by the lattice symmetry. At low tem¬ 
peratures we observe slitting of the ESR absorption line 
due to lifting of the triplet state degeneracy, which is also 
possibly due to the same DM interaction. Additionally 
we observe well resolved ESR absorption lines from the 
inequivalent ladders which allowed an upper estimate of 
interladder exchange interaction. 


II. SAMPLES AND EXPERIMENTAL DETAILS. 

Single crystals of non-deuterated DIMPY were grown 
from the solution by slow diffusion in a temperature gra¬ 
dient. Samples quality was checked by X-ray diffraction 
and magnetization measurements. Concentration of the 
paramagnetic defects estimated from the 500 mK mag¬ 
netization curve is below 0.05%. 

DIMPY belongs to the monoclinic space group 
P2(l)/n with lattice parameters a = 7.504 A, b = 
31.613 A, c = 8.206 A and the angle /? = 98.972° 4 As- 
grown crystals have a well developed plane orthogonal to 
the 5-axis and are elongated along the a direction. 

ESR experiments were performed using set of the 
home-made transmission-type ESR spectrometers at the 
frequencies 18-38 GHz. Lowest available temperature of 
400 mK was obtained by He-3 vapours pumping cryo¬ 
stat. At the measurements below 77 K magnetic field was 
created by compact superconducting magnets, typical 
nonuniformity of the magnetic field at the resonance con¬ 
ditions in our experiments is estimated as < 5... 20 Oe 
depending on the magnet used. High-temperature exper¬ 
iments where done with a resistive water-cooled magnet 
with the held nonuniformity about 5 Oe. 


III. LATTICE SYMMETRY AND POSSIBLE 
ANISOTROPIC INTERACTIONS. 

Monoclinic unit cell of DIMPY includes four magnetic 
Cu 2 + ions that belongs to two spin ladders: two pairs 
of copper ions form rungs of the spin-ladders, which are 
then reproduced by translations along the a-axis. This 
results in the formation of two ladders differently oriented 
with respect to the crystal^ (see Figure [l]). 

Space symmetry of the DIMPY lattice includes inver¬ 
sion center in the middle of each rung and a second order 



FIG. 1. (color online) Crystallographic structure of DIMPY 
with two magnetically nonequivalent spin ladders. Only Cu 
and Br ions are shown along with the main exchange bonds 
Jieg and Jmng■ Solid arrows (red) indicate directions of 
g-tensor main axes for inequivalent ladders, dashed arrows 
(blue) indicate directions of the Dzyaloshinskii-Moriya vec¬ 
tors for inequivalent ladders, as found from the data fit (see 
text). Broad double-headed arrows links DM vector and g- 
tensor axis corresponding to the same ladder. 


screw axis parallel to the crystallographic 5 direction that 
links different ladders. 

These symmetries place strong restrictions on the 
possible microscopic anisotropic interactions in DIMPY 
despite the low crystallographic symmetry. First, all 
anisotropic interactions along the legs of the ladders 
should be uniform because of translational symmetry. 
Second, Dzyaloshinskii-Moria antisymmetric interaction 
Hdm = D • [Si x S 2 ] on the rungs is forbidden by the 
inversion symmetry. The same inversion symmetry re¬ 
quires that direction of the Dzyaloshinskii-Moria vector 
D have to be exactly opposite on the legs of the same 
ladder. Third, inversion center on the rungs of the lad¬ 
der ensures that g-tensor is always the same for the given 
spin ladder so there are no complications of anisotropic 
Zeeman splitting. 

Second order axis establish relations between the g- 
tensor components and Dzyaloshinskii-Moria vector di¬ 
rection in different ladders. In particular, because of 
this second order axis the effective ^-factor values are 
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the same for both ladders for the field applied parallel or 
orthogonally to this axis. 

This analysis neglects anisotropic interladder couplings 
and possible symmetric anisotropic exchange (SAE) cou¬ 
pling Usae = 'hirSiSl. where J MT are compo¬ 

nents of a symmetric exchange tensor A, which is usually 
constrained by condition TrJ MiT = 0. Symmetric inter¬ 
action is allowed both on rungs and legs of the ladder 
and is also constrained by symmetry operations. How¬ 
ever, as we will demonstrate below, our observations 
point that Dzyaloshinskii-Moria interaction is dominat¬ 
ing anisotropic interaction in the case of DIMPY. 

As for the anisotropic interladder couplings, there is a 
possibility that anisotropic couplings between the equiv¬ 
alent ladders stacked in c direction could be important as 
well: Cu-Cu distance in this direction is even less then the 
distance on the rungs (8.2 A against 8.9 A) and suppres¬ 
sion of the Heisenberg exchange interaction in this direc¬ 
tion is most likely due to unfavorable mutual orientation 
of the electron orbitals of bromine ions mediating this 
superexchange route which could be less important for 
the anisotropic spin-spin interactions arising through in¬ 
volvement of differently oriented excited electron orbitals 
mixed with the ground state by spin-orbital interaction^ 
In the present work we neglect this possibility. 

Thus, the main anisotropic interactions in DIMPY are 
really simple to analyse. They include anisotropic g- 
tensor, which is the same for all magnetic ions of the 
given spin-ladder and a Dzyaloshinskii-Moria interaction, 
which is uniform along the leg of the ladder and the 
Dzyaloshinskii-Moria vectors are exactly opposite on the 
legs of the given ladder. 


IV. EXPERIMENTAL RESULTS. 

A. Angular dependence of the ESR absorption at 
77 K. 

We have taken rotational patterns of ESR absorption 
for the magnetic field applied in different crystallographic 
planes. Because of monoclinic lattice symmetry care 
should be taken with consistent determination of the field 
direction with respect to the lattice axes. We will use a 
cartesian basis with X\\a, Y\\b and Z\\c* for the direc¬ 
tion description. Rotation patterns were taken for the 
field confined to {XY) and (XZ) planes and to the plane 
containing Z axis and an (Y — X) direction. All rotation 
patterns were taken for more then 180° angular sweeps. 

Examples of absorption spectra and angular depen¬ 
dences of the ^-factor are shown of the Figure [2] We 
observe one or two Lorentzian absorption lines. These 
absorption lines are clearly due to the different spin lad¬ 
ders: the ladders are equivalent with respect to the mag¬ 
netic field for H||Y and H _L Y, and we observe single 
component absorption at these orientations. Anisotropy 
of (/-factor is typical for the Cu 2+ ion, (/-factor varies 



FIG. 2. (color online) Main panel: angular dependence of the 
(/-factor at 77 K. Symbols - experimental data, curves - uni¬ 
axial (/-tensor best fit (see text). Experimental error is about 
of 0.1% and is within symbol size. Inset: example of ESR 
absorption spectra at representative orientations. Symbols - 
experimental data, curves - best fit with Lorentzian lineshape, 
vertical lines corresponds to certain (/-factor values. Narrow 
line with g = 2.00 (at 6.13 kOe) is a DPPH marker. 


from about 2.03 to 2.30 in agreements with powder ESR 
measurements of Refill 

Angular dependence of the ESR linewidth was deter¬ 
mined by fitting observed absorption spectra with a sin¬ 
gle lorentzian line or with a sum of two lorentzian lines 
(Fig El). Typical half-width at half-height at 77 K was 
around 50 Oe. Field inhomogeneity in the used mag¬ 
net and uncertainties of the fit procedure limit accuracy 
of the linewidth determination to about 5 Oe, however 
angular dependence is clearly present. We were able to 
cross-check our results at certain selected orientations on 
a commercial Bruker X-band spectrometer and we have 
found that X-band data are in agreement with our re¬ 
sults. 


B. Low-temperature ESR. 

Low temperature (below 77 K) ESR absorption was 
measured at certain fixed field directions: for the field 
applied parallel to the symmetry axis H||Y and for the 
field canted by approximately 45° towards X axis. In 
the first case both ladders are equivalently oriented with 
respect to the magnetic field, while the later case corre¬ 
sponds to the maximal difference of the ladders’ effective 
(/-factor, as evidenced by 77 K measurements. 

As expected, we observe single-component ESR ab¬ 
sorption for H||Y and two resolved ESR signals from 
different ladders for the canted sample. Temperature 
evolution of the ESR absorption spectra is qualitatively 
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FIG. 3. (color online) Angular dependence of the ESR 
linewidth at f=17.2 GHz, T=77 K (half-width at half-height). 
Vertical bar at the left panel shows typical errorbar size (dou¬ 
ble error). Squares: width of the high-c; component. Circles: 
width of the low-p component. Curves: model description 
(see text), solid lines show full linewidth, dashed and dot¬ 
ted lines show contributions due to DM and SAE interactions 
respectively. Marks ”1” and ”2” indicate contributions corre¬ 
sponding to the same ladder. 
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FIG. 5. (color online) ESR absorption spectra at the tempera¬ 
ture T « 1 K at different frequencies, H||(A + V). All spectra 
are shifted along the field axis to fit positions of the main ab¬ 
sorption subcomponents. Left panel: left absorption compo¬ 
nent (g « 2.28), weak absorption subcomponent is magnified 
by the factor of 5 or 10 for better presentation. Right panel: 
right absorption component (g « 2.05). Vertical dashed lines 
mark positions of the absorption subcomponent at lowest fre¬ 
quency. Triangles on the right panel mark position of the 
DPPH marker absorption (g = 2.00). 



FIG. 4. (color online) ESR absorption spectra at low tem¬ 
peratures, H||(A + Y). Vertical dashed lines mark resonance 
fields corresponding to the shown (/-factor values. Horizontal 
dashed line at 0.45 K curve is a guide to the eye at zero- 
absorption level. Narrow absorption line at g = 2.00 is a 
DPPH marker. 


similar in both cases (Figure 2]). Below 10 K the ESR 
absorption intensity freeze down due to the presence of 
the energy gap. ESR signal continue to loose intensity 
down to 450 mK and almost vanishes at this temperature. 
Lowest temperature (450 mK) ESR absorption includes 
broad powder-like absorption spectrum probably related 
to the distorted surface of the sample. 


FIG. 6. (color online) Left panel: temperature dependence of 
the ESR intensity below 1 K at f=34.6 GHz, H||(A' + Y). 
Inset: examples of ESR absorption and ESR components 
and subcomponents notations. Symbols - experimental data, 
dashed lines - fits with thermoactivation law I oc exp(—A /T). 
Right panel: dependence of the determined activation gaps for 
different spectral subcomponents. Filled symbols: intense A1 
and B1 subcomponents, open symbols - weak A2 and B2 sub¬ 
components. Lines: parameters-free model dependence cal¬ 
culated with the zero-field gap value A/jvs known from the 
inelastic neutron scattering experiments.—^ Inset: scheme of 
the energy levels of a spin-gap magnet in a magnetic field. 
Solid vertical arrows show transitions corresponding to the 
observed ESR absorption, dashed vertical arrows mark acti¬ 
vation gaps for these transitions. 
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We did not observed any additional absorption signals 
which could be related to the formation of the field in¬ 
duced ordered phase above the critical field to appear at 
the lowest temperature of 450 ml( in the fields up to 10 T 
at the frequencies of 26...35 GHz. This is in agreement 
with the known phase diagram of DIMPY— demonstrat¬ 
ing that highest temperature of the transition into the 
ordered state is about 300 mK. 

Additional splitting of the ESR absorption lines was 
observed around 1 K (Figure [5]), resonance fields of the 
split sub-components differ by approximately 150 Oe. 
This splitting was observed at various frequencies, it 
was most pronounced on the high-held component of the 
canted sample ESR absorption spectra. One of the split 
sub-components is much weaker then the other and freeze 
out faster on cooling. Remarkably, mutual orientation of 
the weaker and stronger sub-components is different for 
the low-held and high-held components. We did not ob¬ 
serve resolved splitting for the H||Y orientation, instead 
a weak peak of the lincwidth was observed around the 
same temperature of 1 K probably indicating unresolved 
splitting. 

At low temperatures intensities of all components fol¬ 
low exponential law I oc exp (—A /T) (Figure [HI). Energy 
gap for the weaker sub-components is larger then that for 
the main subcomponents. By taking temperature depen¬ 
dences of the ESR absorption at different frequencies we 
were able to determine the values of the energy gaps at 
several frequencies revealing dependence of the activation 
energy from the resonance held (Figure [6]). 


C. ESR linewidth evolution from 300 K to 400 mK. 

Temperature evolution of the ESR linewidth was mea¬ 
sured from room temperature down to 400 mK (Figure 
Q. The temperature dependence is qualitatively sim¬ 
ilar in all orientations and demonstrate strongly non- 
monotonous behavior. At high temperatures (above 
90 K) linewidth strongly increases with heating rising 
from about 50 Oe at 77 K to about 300 Oe at 300 K. On 
cooling below 77 K linewidth again increases reaching 
maximum at temperature T max = 9.0 ± 0.2 K. Temper¬ 
ature of the maximum is the same for all orientations, 
while linewidth value at the maximum varies from 90 Oe 
to 140 Oe, both of the extreme values being observed 
in the orientation of maximal splitting H||(X ± Y) for 
different components of ESR absorption. Below T max 
linewidth again decreases reaching its minimal value of 
about 10 Oe (observed at H||Y) at 2 K, which is most 
likely limited by the held inhomogeneity in our setup. 
On cooling below 2 K a peak in the linewidth is observed 
around 1 K. The peak is most pronounced for the high- 
held component in the orientation of maximal splitting 
of the ESR absorption components, peak position coin¬ 
cides with the temperature of subcomponents appear¬ 
ance. Similar but less pronounced peak is observed for 
H||Y. High-g component in the H||(X + Y) orientation 
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FIG. 7. (color online) Temperature dependence of the ESR 
linewidth. Upper panel: H||(A' + Y). Circles: low-held (high- 
g) component, squares: high held (low-//) component. All 
data are measured on the samples from the same batch. Lower 
panel: H||Y. Filled and open symbols corresponds to the data 
measured on samples from different batches. Experimental 
data were collected in different experimental setups operat¬ 
ing at different frequencies, vertical lines mark approximate 
temperature boundaries for different experiments, microwave 
frequencies of each experiment are given. Typical errorbar 
size is around symbol size. Curves on both panels show em¬ 
pirical ht equations (see text and Table [Hi. 

does not demonstrate such a peak, which is probably re¬ 
lated to the very low intensity of the appearing weaker 
subcomponent which cause fitting procedure to lock on 
the main spectral subcomponent. Finally, on cooling be¬ 
low 1 K linewidth of both components in the H| | {X + Y) 
orientation increases again. 


V. DISCUSSION. 

A. Recovery of the (/-tensor. 

Observed angular dependences of the (/-factor can be 
fitted assuming uniaxial (/-tensor. As was described in 
the Section uni (/-tensor is the same for the given lad¬ 
der and orientations of the (/-tensors in inequivalent lad¬ 
ders are bound by 2-nd order axis. Hence, directions 
of the main axis can be expressed via polar angles as 
n g i , 2 = (± sin O cos (/>; cos 0; ± sin 0 sin <f>), here we count 
polar angle 0 from the 2-nd order axis Y||6, different 
signs corresponds to the different ladders. 

Least squares fit of our data (see Figure 0) yields g- 
tensor components g\\ = 2.296 ± 0.010 and g± = 2.040 ± 
0.006 and angles 0 = (34.8 ± 1.5)° and (/> = (178 ± 4)°. 
Fit quality can be improved by assuming general form of 
the (/-tensor. However, this results only in minor planar 
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anisotropy (with principal ^-factor values of 2.038±0.010 
and 2.058 ± 0.010) which is on the edge of experimental 
error. 

Main axis of the g-tensor within accuracy of our exper¬ 
iment lies in the (XY)-plane of the crystal. This seems to 
be accidental as there is no symmetry reasons to choose 
this plane in the monoclinic crystal. Orientation of the 
^-tensor main axes with respect to the crystal struc¬ 
ture is shown at the Figure[T] We can not decipher which 
of the orientations corresponds to different ladders. 

Found values of the main g-tensor components coincide 
with the values found in earlier powder high-field ESR 
experiment^ The value of the g -factor for the H||a case 
g a = 2.17 was found in Refflll by magnetization fit, this 
value disagree by 2% with the value g = 2.130 ± 
0.005 found in our experiment. 


B. Inter ladder coupling estimation. 

Anisotropy of the <?-tensor opens a direct way to the es¬ 
timation of interladder coupling. If the coupling between 
the ladders with different g-tensor orientation would be 
strong enough, then a common spin precession mode 
would be observed, a well known exchange narrowing 
phenomenom 15 i 16 

Instead, we observe well resolved ESR absorption sig¬ 
nals from the inequivalent ladders. Thus, upper limit on 
the interladder coupling can be estimated from the min¬ 
imal splitting AH observed, which is around 40 Oe at 
17 GHz experiment (corresponds to components g -factor 
difference A g = 0.015 at Figure 0). 

As it is known from the exchange narrowing theory, 
equally intense components with splitting Aw will form 
a common precession mode if the coupling strength is 
above J c ~ frAui. Hence we obtain an estimate Jinter < 
~ 5 mK or about 0.5 geV. 

This value is in reasonable agreement with the ear¬ 
lier estimate of the interladder exchange coupling from 
the ordered phase boundary calculated in mean-field 
approximation^ as nJ' MF = 6.3 peV (here n is the num¬ 
ber of coupled ladders, coupling considered to be equal in 
all directions). Note that our observation provides a di¬ 
rect estimate of the coupling between the unequivalently 
oriented ladders only. 


C. Linewidth temperature dependence. 

Non-monotonous temperature dependence of the 
linewidth indicates that spin precession relaxation is gov¬ 
erned by different processes in the different parts of the 
studied temperature range. We fit this temperature de¬ 
pendence by set of empirical equations as shown on Fig¬ 
ure [7] and discussed below. Values of the fit coefficients 
for the empiric equations used are gathered in the Table 

m 


High temperature increase of the ESR linewidth (above 
77 K) is naturally related to the spin-lattice relaxation: 
increase of the phonons population numbers leads to the 
increase of the relaxation rate. Linewidth dependences 
can be fitted by the sum of the constant contribution 
describing the high-temperature spin-spin relaxation and 
an empirical activation law AH = AHo+Aexp(—E a /T). 
Activation energy is E a = (1400 ± 150)X. Similar be¬ 
havior with activation energy of the same order of mag¬ 
nitude was reported for other cuprate s 17 ! 18 and it was 
discussed^ as a relaxation via excited state with a com¬ 
peting Jan-Teller distortion. However, detailed analy¬ 
sis of the lattice relaxation is beyond the scopes of the 
present paper. 

Lattice contribution vanishes with cooling. Shallow 
minimum of the linewidth at 70-100 K indicates that 
phonon relaxation channel is practically frozen down 
here. Hence, we assume that linewidth measured at 77 K 
is mostly due to spin-spin relaxation. 

It is well know n 19 : 20 that anisotropic spin-spin in¬ 
teractions are responsible for the spin-spin relaxation. 
Thus, ESR linewidth provides access to determine these 
interactions strength. For the concentrated magnets 
Dzyaloshinskii-Moria interaction (which is allowed by the 
lattice symmetry of DIMPY) and symmetric anisotropic 
exchange interaction are the main contributions. Tem¬ 
perature dependence of the ESR linewidth is one of the 
physical effects to find which of theses anisotropic inter¬ 
actions dominates the linewidth. 

Detailed description of the ESR linewidth in a quan¬ 
tum spin-ladder is only emerging now: case of spin lad¬ 
der with symmetric anisotropic spin-spin coupling was 
considered by recently by Furuya and Sato^i, a theory 
accounting for the uniform Dzyaloshinskii-Moriya inter¬ 
action is still to be constructed. Theory of an ESR 
linewidth for a quantum S = 1/2 chain was devel¬ 
oped by Oshikawa and Affleck more then decade a,goM 
We will apply their results to understand qualitatively 
temperature dependence of the ESR linewidth at high 
temperatures T Ji eg , r ung■ Oshikawa and Affleck 
have demonstrated that contribution of the symmetric 
anisotropic exchange interaction to the ESR linewidth 
(exchange anisotropy in their terms) decreases with cool¬ 
ing. They also considered contribution of the stag¬ 
gered Dzyaloshinskii-Moria interaction and have found 
that in this case ESR linewidth is increasing as 1/T 2 at 
low temperatures. Oshikawa and Affleck demonstrated 
that at high temperature limit staggered Dzyaloshinskii- 
Moria interaction results in the linewidth increasing with 
cooling as 1/T. As the high-temperature linewidth is 
determined by the pair spin correlations, this result 
should be actually the same for the staggered and uni¬ 
form Dzyaloshinskii-Moria interaction. This conclusion 
is in agreement with the results of Refill for the uni¬ 
form Dzyaloshinskii-Moria interaction in quasi one di¬ 
mensional antiferromagnet CS 2 CUCI 4 . Thus, increase of 
the linewidth with cooling below 80 K is a direct indica¬ 
tion of the dominating role of the Dzyaloshinskii-Moria 
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TABLE I. ESR linewidth empirical fit equations parameters in different temperature ranges and at different orientations. For 
the orientations with two ESR components resolved (LF) and (HF) marks fit results for the low-field and high-field components 
correspondingly, (HF+LF) marks cases where linewidths of both components are close and their temperature dependences are 
fitted jointly. 


Temperature range 
and fit eqns. 

Parameters 

Typical fit 

accuracy 

H||Y 

H||(x ± y) 

H||(f/ ± z) 

T > 80 K 

AH 0 , Oe 

±5 Oe 

46 

45 (HF); 56 (LF) 


AH = AHo +Aexp{-E a /T) 

A, xlO 4 Oe 

±50% 

5.9 

1.7 (HF); 1.8 (LF) 



Ea, K 

±150 K 

1510 

1240 (HF); 1310 (LF) 


15K< T < 80 K 

Affoo, Oe 

±5 Oe 

35 

35 (HF); 44(LF) 

40 (HF±LF) 

AH = Ai7oo(l + 0/T) 

0, K 

±3 K 

24 

16 (HF); 22(LF) 

14 (HF±LF) 

1.5K< T < 7 K 

AH 0 , Oe 

±2 Oe 

10 

15 (HF);15 (LF) 

15(HF±LF) 

AH = AHo +Aexp(-E a /T) 

A, xlO 2 Oe 

±30% 

6.4 

1.7 (HF);3.9 (LF) 

5.3 (HF±LF) 


Ea, K 

±1.5 

14 (9.6^) 

6.4 (HF);8.4 (LF) 

10.5(HF±LF) 


a Value of E a = 9.6 K corresponds to the best fit with fixed parameter All o = 0 (see text). 


interaction for the spin relaxation processes in DIMPY. 
To model first order of the 1 /T expansion we fit our data 
by the law AH = Ai? 00 (l ± 0/T). Characteristic tem¬ 
perature 0 is anisotropic and varies from 15 to 25 K in 
the orientations presented on the Figure [7] This tem¬ 
perature scale is close to the exch ange i ntegral value in 
agreement with the results of Refs l22lll7l 

The crude estimate of this interaction strength can 
be obtained from the linewidth at 80 K, which is ap¬ 
proximately 50 Oe. As this temperature far exceeds the 
exchange integral scale high-temperature approximation 
can be used. We will discuss exact calculations below 
while describing angular dependence, but as an estimate 
one can write TiAlo ~ EL. or D ~ gjinAHJ ~ 0.3 K. 

As the temperature decreases below approximately 
10 K linewidth start to decrease. This decrease is natu¬ 
rally related to the gapped spectrum of the spin ladder. 
At low temperatures magnetic properties of a spin ladder 
can be described on the triplet quasiparticles language 
and linewidth is then interpreted as an inverse lifetime 
of these quasiparticles, which is partially determined by 
their interaction. As temperature approaches scale of the 
energy gap, quasiparticles population numbers decreases, 
gas of the quasiparticles became diluted and quasiparti¬ 
cles interaction contribution froze out. This results in the 
narrowing of the ESR absorption line with cooling. The 
linewidth temperature dependence indeed follows ther¬ 
moactivation law AH = AHq + Aexp(—E' a /T) with ac¬ 
tivation energy E' a = 6.4... 14 K in different orientations. 

If the relaxation processes would be due to the pair in¬ 
teraction of quasiparticles, the relaxation rate would be 
proportional to the quasiparticles concentration squared 
and the activation energy would be about 2Ao, where 
A 0 = 0.33 meV (or 3.8 K) is a zero-field gap. The ac¬ 
tivation energies for H||(X ± Y) orientation are close to 
this expectation. However, for the field applied along 
the second-order axis H||Y the best fit activation energy 
(14 ± 1.5 K) far exceeds the doubled zero-field gap value. 


Most likely this result is an artefact due to the effects of 
field inhomogeneity in our experimental setup: the low 
temperature linewidth is minimal for H||Y and could be 
limited by experimental resolution. This leads to over¬ 
estimation of AHq parameter which in turn results in 
overestimation of the activation energy. Tentative fit of 
the H||Y data with AHq value fixed to zero yields ac¬ 
tivation energy of 9.6 K which is much closer to twice 
zero-field gap value. Similar AHq = 0 fits in other orien¬ 
tations lead to the smaller corrections of the determined 
activation energy. 

Thus, within accuracy of our experiment, which is 
mostly limited by field inhomogeneity of the magnetic 
field in the superconducting coil used, we can con¬ 
clude that pair quasiparticles interactions dominates 
spin-spin relaxation processes of the spin-ladder in low- 
temperature regime. 

The peak of the linewidth around 1 K is related 
to the splitting of the ESR lines into subcomponent, 
its origin is related to the classical exchange narrow¬ 
ing phenomenom 15 i 16 The exchange frequency became 
temperature dependent being related to the quasiparti¬ 
cles concentration. At low temperatures (low quasiparti¬ 
cles concentration) split ESR line is observed, at higher 
temperatures (higher quasiparticles concentration) effec¬ 
tive exchange interaction between the quasiparticles gain 
efficiency and a common precession mode is formed. 
Crossover between these regimes results in the broad¬ 
ening of ESR line. Similar effect is observed in other 
spin-gap magnet o 23 : 24 and in other systems^ 

Finally, definitive increase of the linewidth below 
700 mK is probably indicative of the critical regime in 
the vicinity of the field induced phase transition. 












D. Angular dependence of the ESR linewidth. 

According to the theory of the exchange narrowed res¬ 
onance spectr a 26 ! 27 , the half width at half maximum for 
a single Lorentzian shaped line is given by 


Here angles a and /3 define orientation of the local co¬ 
ordinate system {x, y, z) where Zeeman term in Eq. ([1]) 
takes diagonal form g/igHS 2 with 

g = VA 2 + B 2 ±C 2 , ( 8 ) 


AH = C 


\ Mi 


1/2 


M 4 


( 2 ) 


where C is dimensionless constant of order unity, depend¬ 
ing on how the wings of Lorentzian profile drop at fields 
of the order of exchange field {J/gg b "C A7L)21; M 2 and 
M 4 are the second and fourth moments of resonance line, 
firstly introduced by Van Vlecb^ 


M 2 


<[< H anis ,S+}{S-,H anis ]) 


h 2 {S+S~) 


( 3 ) 


M 4 = 


([Hex, [H anis , 5+]] [[S~ , Hants] ,U B x\) 

hA(S+S~) 


( 4 ) 


where S ± denote left/right circular components of the 
total spin summed up over the whole sample, H ex is 
isotropic exchange Hamiltonian, Hants is anisotropic one 
that doesn’t commute with H ex hence causing broaden¬ 
ing of the resonance line. 

Analysis of the ESR linewidth based on calculation of 
the spectral moments is a well developed method which 
allows to identify nature of spin-spin interactions and 
estimate their magnitudes in magnetically concentrated 
systems^ Its benefit is that in high temperature limit 
(T -A oo) an exact expression for linewidth can be found 
out for an arbitrary spin system, whatever space dimen¬ 
sion and exchange couplings ! 20 ! 30 

In the present paper we apply the ’’method of mo¬ 
ments” to a strong-leg spin ladder system, described by 
Hamiltonian © with uniform Dzyaloshinskii-Moriya in¬ 
teraction 


Hdm = EE Dz[S* 5 j x Sz,i+i], (5) 

i 1= 1,2 


here i enumerates rungs of the ladder and l enumerates 
legs of the ladder, DM vectors on the legs are considered 
arbitrary for the moment (D 4 =/= D 2 ). Substituting Eq. 
(ED into Eqs. ED, ED and using the corresponding com¬ 
mutation relations for S' = 1/2 spin operators, for the 
linewidth Eq. ED i n high temperature limit we have 


A Hg M {Oe) = C £ 

1 = 1,2 


[Pg±Dg±2P 2 ] ; 
iV^n B JDg(0, <j>) 


( 6 ) 


where Jd = Jf e g + 2 J? U ng mean an average exchange 

integral , 30 and angular dependence is determined by the 
transformation for DM vector 

D x = Dx cos f3 cos a + Dy cos /3 sin a — Dz sin f3, 

D y = Dy cos a — Dx sin a, (7) 

D z = Dx sin /3 cos a + Dy sin /3 sin a + Dz cos /3. 


where 

A = g X x sin 0 cos <f> + g X y cos © + g X z sin 0 sin</>, 

B = gyx sin 0 cos 4> + <Jyy cos 0 + gyz sin 0 sin <j>, 

C = gzx sin 0 cos + gzy cos 0 + gzz sin 0 sin <j> 
and 

A 0 C 

cos a= . =, cos p = - . = = =, 

VA 2 + B 2 VA 2 +B 2 + C 2 

here polar 0 and azimuthal <f> angles define the direction 
of external magnetic field, so that 0 and </> are counted 
from Y and X axes, respectively, as during the g-tensor 
recovery procedure. 

Note that by setting J rung = 0 and D 4 = D 2 in Eqn. 
© we immediately arrive to the known result for ID 
Heisenberg chain with uniform DM interaction [see for¬ 
mula (16) in Ref. Qj}. 

DIMPY has two inequivalent ladders with different g- 
tensors and DM vectors. In accordance with crystal sym¬ 
metry of the DIMPY (see Sec. Hill) . the legs within same 
ladder are linked by inversion, so that 


g( fc) = g { 2 k \ D( fe) = -D^\ (k= 1 or 2), (9) 

here upper index (k = 1 , 2 ) denotes the inequivalent lad¬ 
ders and lower index enumerate legs od the ladder. The 
legs of inequivalent ladders are linked by screw rotation 
along the second order axis, hence 

d[ 2) =c , 2 (y)d| 1) , 

* ( 10 ) 

gp ) =C 2 (Y)g i ( 1 ) C' 2 (Y)- 1 , (l = 1,2), 


Orientation of the local axes is essentially different for 
inequivalent spin ladders. Having known directions of 
the main axes of 3 -tensors (see Sec. IV A[) . it’s easy to 
find their components referred to crystallographic axes 


5 (1,2) = 


1 2.128 =f0.12 -0.008' 

=F0.12 2.214 ±0.005 

y—0.008 ±0.005 2.038 


( 11 ) 


where upper(down) sign corresponds to the ladder with 
upper (down) sign of n g . 

Simulation of experimental data on the linewidth 
angular dependence by Eq. © showed that the 
Dzyaloshinskii-Moria interaction describes the angu¬ 
lar variation of the linewidth in DIMPY well enough 
(within experimental error). However model including 
Dzyaloshinskii-Moria interaction only predicts value of 
linewidth which is systematically less then the exper¬ 
imental values by about 12 Oe. This fact indicates 
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there is an additional (small compared to DM inter¬ 
action) source of the line broadening in DIMPY. The 
modelled values of the linewidth can be reconciled with 
the experimental ones by adding isotopic contribution 
AHq = 12 Oe, which can be probably ascribed to the 
residual spin-lattice relaxation, or by considering other 
anisotropic spin-spin couplings. Contribution to the 
linewidth from dipole-dipole interaction is quite small for 
DIMPY and at the shortest distance between Cu ions 
(r = a « 7.5 A) following conventional estimation^ it 
does not exceed ~ 0.5 Oe. Additional broadening in 
DIMPY can be related to SAE interaction along legs and 
rungs of the spin ladders which usually appear as further 
sources of ESR line broadening beyond the dominant DM 
interactional U3£ The contribution to the linewidth due to 
SAE interaction is derived in Appendix lAl 

Taking into account symmetry relations (Eqs. (fUjl, 
(HUD, (ED, (1A3I) ~1 and, for definiteness setting D ^ 1 1 = D, 
= A, the fitting of linewidth angular dependence 
yields D x = 0.15, D v = -0.14, D z = 0.075 K and al¬ 
most diagonal exchange-tensor with components Jxx = 
0.08, Jyy = —0.03, Jzz = —0.05, Jxy = —0.015 K and 
Jxz = Jyz = 0. During simulation the Lorentzian pro¬ 
file with exponential wings was assumed, which implies 
C' = 7t v / 2. 

Directions of the found D vectors are shown on 
the Figure |T| Components of the DM vector 
and SAE tensor given above correspond to the lad¬ 
der with the g-tensor main axis orientation n g i = 
(sin0cos<)>;cos0;sin0sin(/>) r; (—0.57; 0.82; 0) (when 
comparing with FigQ] note that n g and — n g are phys¬ 
ically equivalent), angle between D and n g vectors (re¬ 
duced to (—f; f) range for convenience) is approximately 
23°. Magnitude of the obtained Dzyaloshinskii-Moria 
vector |D| = 0.22 K agrees well with the crude estimation 
above (section IVCl) . 

As it is seen from Fig. [UJ taking into consideration 
only an exchange mechanisms of spin anisotropy within 
the legs of ladders gives a good compliance with experi¬ 
ment. However, it is necessary to stress out, that without 
DM interaction SAE coupling only (see Appendix [Aj to¬ 
tally failed to give a correct description of the angular 
dependence of linewidth in DIMPY. 

Our simulation shows that absolute value as well as 
angular anisotropy of the linewidth are predominantly 
attributed by DM interaction, while contribution to the 
linewidth due to SAE interaction is relatively small com¬ 
pared to DM one and gives a weak angular variation. 
Similar behavior of ESR linewidth with coexistent con¬ 
tributions from DM and SAE interactions within 5 = 1/2 
antiferromagnetic chains was observed in hight symmetry 
crystal structure KCuF^. 

The found DM vector have not only transverse but 
also nonzero longitudinal (with respect to the Cu-Cu 
exchange bond) component within the legs. Such re¬ 
sult does’t contradict with the general rules for DM vec¬ 
tor, established by Moriya-^ based on general symmetry 
grounds for a pair of exchange interacting ions. More¬ 


over, a simple analysis of the recovered g-tensors (see Sec. 
IV All leads to the same conclusion about direction of the 
DM vector. Since the axial component of the g-tensor 
has a maximal value, then the ground state orbitals of 
Cu 2+ ions (typically ”i 2 — g 2 ”-like symmetry and 5||n s ) 
should predominantly lie within the plane perpendicu¬ 
lar to the main axis of a g-tensor, because the maximal 
matrix element (< x 2 — y 2 \l g \xy >= —2i) relevant to 
spin-orbital coupling appears only in the case when an 
external magnetic field is applied parallel to the main 
axis of g-tensor. For the same reason an effective DM 
vector predominantly should lie along the main axis of 
the g-tensor. It should be noted that conventional rule 
determining DM vector as D oc [ni x 772 ] i 34 ’ 35 where n\ 
and 712 are unit vectors connecting a exchange interacting 
ions with bridging ion, is not applicable in present case. 
Possible failure of this rule was mentioned before in Ref. 
[36l . referring to the features of exchange process through 
a two bridging ions, which is also the case of DIMPY (see 
Fig. [ID. 

Thus, analysis of ESR linewidth allowed us to conclude 
that the DIMPY is a rare case of compound in which DM 
vector has a component along the line connecting the pair 
of exchange interacting ions. This is a consequence of 
low crystal symmetry of DIMPY and nontrivial orbital 
ordering. 


E. Low-temperature sub-components appearance. 

First, we recall main observations on the subcompo¬ 
nent appearance. ESR components splits around 1 K 
into two sub-components, one of which is much weaker. 
The splitting is best observed at H||(X + Y) orientation. 
Position of the weaker sub-component with respect to 
the stronger sub-component is different for both ESR ab¬ 
sorption components. Maximal splitting is about 150 Oe 
and it decreases as the resonance field approaches critical 
field, weaker subcomponent became unresolvable at the 
fields above 2/3 of the critical field. Activation energies 
for the stronger and weaker subcomponents are different. 

All these observations can be explained as an effect 
of the zero field splitting of triplet sublevels. This ef¬ 
fect was already observed for various spin-gap magnets, 
e.g. TlCuCl3f22 or PHCCi^i Anisotropic interactions lift 
degeneracy of the 5 = 1 triplet state and frequencies 
of the dipolar transitions | S z = +1) «-)■ \S Z = 0) and 
|5 Z = — 1) «-)■ |5^ = 0) would become different. Here we 
assume, which is perfectly valid for the case of DIMPY, 
that the anisotropy is very small and spin projection 
on the field direction S z can be considered as a good 
quantum number. Therefore, in the presence of such an 
anisotropy the resonance fields for |5 Z = +1) -O- |5 2 = 0) 
and |5 Z = —1) O \S Z = 0) transitions in the constant 
frequency ESR experiment would differ and ESR absorp¬ 
tion split into two sub-components. 

Observed difference of the activation energies for the 
absorption sub-components and dependence of the acti- 
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vation energy on the microwave frequency used in the 
experiment is a direct consequence of this explanation. 
The ESR intensity at low temperature is determined by 
the population of the lowest sublevel. Hence, for the 
| <S Z = +1} O |S Z = 0) transition the activation en¬ 
ergy is A ss Ao — gnBH res = Ao — his, being deter¬ 
mined by the population of the |S* = +1) sublevel (en¬ 
ergy of this sublevel decrease with field, see inset on Fig¬ 
ure 0. In the same time the activation energy for the 
|S* = — 1) <->■ |S* = 0) remains constant (and equal to 
Ao) since the energy of |S* = 0) sublevel is field inde¬ 
pendent. The dependences of the activation energy on 
the microwave frequency of the ESR experiment are de¬ 
scribed by this model parameter-free using the zero-field 
gap value of 0.33 meV from the inelastic neutron scatter¬ 
ing experimented^ 

Behaviour of the sublevels of the spin-gap magnet in 
the vicinity of the critical field is a long-discussed prob¬ 
lem. There is a general macroscopic (or bosonic) ap¬ 
proach of Refs. U [Bj and a ID fermionic approach of 
Tsveliki developed for the spin-chains. Fermionic model 
of Tsvelik yields results formally equivalent to the results 
of perturbation treatment of anisotropic interactions 
Thus, within these approaches the sublevels behave lin¬ 
early in the vicinity of the critical field and the splitting of 
the ESR subcomponent should be then field independent. 
Bosonic model, on the contrary, predicts square-root-like 
approach to the critical field for the low-energy sublevel, 
while field dependence of the high-energy sublevel re¬ 
mains linear in the vicinity of the critical field. Therefore, 
sub-components splitting will change close to the critical 
field. However, this nonlinearity of the bosonic model 
extends only in the small vicinity of the critical held 
( H c — H) ~ AE/hb ~ A H, here A E is the zero held 
triplet sublevels splitting and AH ~ 150 Oe is the ob¬ 
served sub-components splitting. We observe (Figure [5]) 
that the observed splitting is halved (compare 50.18 GHz 
and 26.30 GHz curves at the Figure 0 in the held of 
about 2/3 of the critical held (H c ~ 30 kOe, zero-held 
gap of 0.33 meV corresponds to the frequency of 80 GHz), 
i.e. well below this nonlinearity range. This probably in¬ 
dicates that held evolution of the split sub-components 
follows some other laws on approaching the critical held. 
Similar behaviour of the ESR line split by the uniform 
Dzyaloshinskii-Moria interaction was recently reported 
for a quasi-ID antiferromagnet Cs 2 CuCl 4 ^ 

Under an assumption that the uniform Dzyaloshinskii- 
Moria interaction along the legs of the ladder is respon¬ 
sible for the observed splitting, anisotropy axis have to 
be aligned along the D vector. We calculated effects of 
the DM coupling perturbatively for the limiting case of 
strong-rung ladder (see Appendix iBl). Interdimer DM in¬ 
teraction mixes one and two-particle excited S z = ±1 
states which results in the triplet sublevels splitting by 
SE = yj, S z = ±1 sublevels being shifted down. This 
corresponds to the easy-axis anisotropy for the triplet 
excitations, D direction being the easy axis direction. 
Taking the magnitude of the DM vector D « 0.20 K 


as estimated from the high-temperature ESR linewidth 
analysis and substituting energy gap of 0.33 meV as an 
exchange parameter of the perturbative model we ob¬ 
tain an estimate of the sublevels splitting 5E ~ 5 mK 
which corresponds to sub-components splitting of about 
40 Oe, factor of four less then the experimentally ob¬ 
served value. However, perturbative treatment starting 
from the uninteracting dimers is at best a qualitative 
model for a strong-leg ladder and a detailed description 
of a strong-leg spin ladder with uniform Dzyaloshinskii- 
Moria interaction needs a separate theoretical effort. 

We can not unambiguously determine type of the 
anisotropy from our experimental observation since our 
setup does not allow to take an angular dependence 
at He-3 temperature range. However, as it is known 
from the formally similar problem of S' = 1 ion in 
a crystal fiel d 19 ’ 20 the effective anisotropy constant 
changes monotonously with field rotating away from the 
anisotropy axis C e f f — % {3 cos 2 £ — l), where £ is an 
angle counted from the anisotropy axis z and anisotropy 
C enters spin Hamiltonian as C(S Z ) 2 . It is maximal at 
the field parallel to the anisotropy axis, it change sign 
and decreases by the factor of two at the orthogonal ori¬ 
entation of the magnetic field and it turns to zero at a 
magic angle. Thus, as splitting observed for the high-held 
component is larger then that from the low-held compo¬ 
nent (approximately 150 Oe vs. 110 Oe, see Figure 0 
and weaker sub-component is located on different side 
from the main sub-component, we hnd it more likely that 
the high-held component corresponds to the ladder with 
the magnetic held close to the true anisotropy axis. In 
this case, as for the held applied close to anisotropy axis 
the weaker subcomponent is located to the right from 
the stronger subcomponent, the splitting of the triplet 
sublevels follows easy-axis type of anisotropy energy of 
S z = ±1 states being lower then energy of S'* = 0 state 
in zero held. 

However, this tentative identification of the anisotropy 
axis deviates from the simple model of DM interaction 
only: as the vectors D and n 3 are quite close for the 
given ladder the low-held component (corresponding to 
the higher longitudinal g-factor) should then be closer to 
the anisotropy axis. Possible reason for this deviation is 
the effect of symmetric anisotropic exchange on the rungs 
of the ladder (see Appendix [B]). S AE coupling is smaller 
in magnitude, but it enters to the triplet sublevels split¬ 
ting linearly, why DM contribution is quadratic. This is 
contrary to the linewidth calculations where both cou¬ 
plings enter quadratically. Thus, description of the sub¬ 
components splitting probably lies beyond the simple 
model with DM interaction only and requires account¬ 
ing for other anisotropic interactions. 


VI. CONCLUSIONS. 

The strong-leg spin ladder system DIMPY is an estab¬ 
lished test example of the Heisenberg spin ladder. How- 
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ever, anisotropic spin-spin interactions, and in particular 
Dzyaloshinskii-Moria interaction of intriguing geometry: 
uniform along the leg of the ladder and exactly oppo¬ 
site on the other leg, give rise to a family of interesting 
phenomena. 

We have estimated parameters of Dzyaloshinskii-Moria 
interaction from high-temperature data. We observe 
splitting of the ESR line at low temperatures which is 
related to the zero-field splitting of the triplet sublevels 
by the same interaction. Finally, we observe series of 
crossovers between different regimes of relaxation of spin 
precession on cooling from room temperature to 400 mK. 

We present qualitative explanations of our observa¬ 
tions. Simple geometry of the exchange couplings and 
anisotropic spin-spin interactions makes DIMPY one of 
the few candidates for the model-free microscopic de¬ 
scription of the effects of anisotropic interactions on the 
properties of a spin-gap magnet, which is still awaiting 
for a theoretical effort. 
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Appendix A: Linewidth contribution by a SAE 
interaction along the legs. 

Symmetric anisotropic exchange (SAE) coupling is al¬ 
lowed both on the rungs and on the legs of the lader. 
Our aim here is to demonstrate that SAE coupling can 
explain contribution of about 20 % of the total linewidth 
that can not be described by DM coupling alone. We 
will focus here on a SAE coupling along the legs of the 
ladder. We have checked that a SAE coupling along the 
rung yields similar angular dependence and its contribu¬ 
tion differs only by some numerical scaling factor. How¬ 
ever we expect that contribution of the SAE couplings 
on the rungs should be small in a strong-leg ladder since 
the overlapping of the orbitals along the rung is smaller. 

An expression for ESR linewidth due to a SAE cou¬ 
pling along the legs of spin ladder is derived similarly as 
was done for a DM one in Sec. IV D1 applying 

TIsae = EE Sz,;AiSz, l+ i (Al) 

i 1= 1,2 

to the Eqs. ©, ®, © in the framework of ’’method of 
moments”, that in the high temperature limit of S = 1/2 
leads to 


AH, 


SAE 


= C 


[(2A ZZ - A xx - X yy ) 2 + 10(A^ Z + A 2 ) + (A xa; - A yy ) 2 + 4A 2 ]i 


1 = 1,2 


8VQ^BJsg(S,4>) 


(A2) 


where Js = J? eg + 2/3 J 2 ung , and the exchange-tensor 

components in the local coordinates, A vl (a,/3), (p, 7 = 
x 1 y , z) are defined by transformation given in Ref. |l3| [see 
formulae A2-A4 therein] and, as usually, are expressed 
via the intrinsic exchange parameters J MT of symmetric 
tensor A in crystallographic coordinates (/x, r = X,Y, Z). 

In accordance with crystal symmetry of the DIMPY 
(see Sec. IIIII) . symmetrical exchange tensors on the legs 
of same ladder are equal, while on the legs of inequivalent 
ladders they are related by a screw rotation along the 
second order axis, that is 

A^=A< fe) , (A = lor 2) 

Aj 2) =G 2 (Y)A' 1) C' 2 (Y)- 1 , (l = 1,2) 

Generally, taking into account relation of J v . y = J rn for 


a symmetric tensor, an exchange tensor has six different 
components, constrained by TrJ = 0 condition. How¬ 
ever, since the anisotropy of the g-tensor and the tensor 
of SAE coupling originates from the same spin-orbital 
interactions, accidental smallness of gxz and gyz com¬ 
ponents allows to assume Jxz = Jyz = 0 during the fit¬ 
ting procedure. We have found, that the remaining four 
components of the SAE tensor are enough to reproduce 
our data. 
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Appendix B: Perturbative treatment of triplet 
sublevels splitting by uniform DM and SAE 
couplings. 

Strong-rung J rung Ji eg limit allows to use wavefunc- 
tions of the isolated dimers as a zero-order approxima¬ 
tion. For the dimer located at the n-th rung wavefunc- 
tion of the ground state is o = ^(| tl)n — | 4-t)n) and 
wavefunctions of the excited triplet are ip n ii = | tt)n, 
'iftnio — ~A 2 (I I 4'T)n) and ip n i—i — | Wave- 

function of the collective ground state is 

and single-particle excited states with the excitation at 

the n-th dimer can be build as 

n— 1 N 

^ill = ^pO x iftnll x ^pO (Bl) 

p—0 p—n+1 

and similarly for other spin projections. Manyparticle ex¬ 
cited state can be constructed similarly keeping in mind 
hard-core repulsion as only one excited state per dimer 
is allowed. 

One-particle states are TV-fold degenerated, this degen¬ 
eration will be lifted by interdimer exchange coupling Ji eg 
giving rise to excitations dispersion. 

We consider effect of interdimer DM interaction m 
with DM vectors oppositely aligned on the legs of the lad¬ 
der. This configuration conserves symmetry axis (direc¬ 
tion of the DM vector, which wc will use as ^-direction), 
thus excitations will have well defined S z values. Inter- 
dimer DM interaction can then be expressed as 

- c+ _ q+ q- _ 
n 

~^2,n^2,n+1 + ‘S'^rA.n+l] (B2) 

By applying this operator to the ground state and to 
one-particle excited states we obtain 

F^ (0) = 0 (B3) 

V^o = 0 (B4) 

1 = |(-< 2 1 ) 0;(„ + 1)11 + *2o;(n-l)ll) ( B5 ) 

( 2 ) 

here dyT/ 0 . rnll are two-particle excited states with S = 
1, S z = 0 excitation on the n-th rung and S = 1, S z = 1 
excitation on the m-th rung. 

Thus, interdimer DM interaction mixes S z = ±1 single 
particle excited states with S z = ±1 two particle excited 
states. This mixing results in the second-order perturba¬ 
tive correction to the energy of the single-particle state 


n 2 

SE = —^rr —, S z = ±1 being shifted down. As this 
correction is the same for all S z = ±1 states, weak inter¬ 
dimer Heisenberg coupling Ji eg will not affect it. 

This result differs from the effect of intradimer DM in¬ 
teraction which mixes £ = 1,5* = 0 state with 5 = 0 
state and shifts the energy of the S = 1 ,S Z = 0 up by 
~r^j ~—. However both interdimer and intradimer DM in- 

rung 

teraction results in the easy-axis anisotropy for the triplet 
excitations (energy of the S z = ±1 states is lower then 
the energy of the S z = 0 state). 

Note also, that the effective anisotropy for the triplet 
excitations is easy axis, while usually anisotropy due to 
DM interaction (e.g., anisotropy of the order parameter 
in the ordered state of an antiferromagnet) is of easy 
plane type. This ’’inversion” of anisotropy seems to be 
a common feature of all spin-gap magnets: it was ob¬ 
tained by perturbative analysis of the role of single-ion 
anisotropy in Haldane magnet^ and was observed in a 
Haldane magnet PbN^A^Os^ similar inversion of the 
anisotropy type between the anisotropy of triplet excita¬ 
tions and order parameter anisotropy in a field-induced 
ordered phase of a spin-gap magnet follows from macro¬ 
scopic approach^. 

Effect of the symmetric anisotropic coupling on a 
strong-rung ladder can be considered similarly. Our aim 
here is to illustrate that its contribution is linear on cou¬ 
pling parameter along the rung so we consider simple 
axial SAE coupling in the form 


F = A leg ^ S ? n £? n+1 + A rung ^ £?,„££,„ (B 6 ) 

n,j n 

here n enumerates rungs (dimers) and j enumerates legs 
of the ladder (spins in the dimer), Ai eg and A rung are 
SAE coupling constants along the leg and rung of the 
ladder. 

By applying this operator to the ground state and to 
the one-particle excited states one can ascertain that SAE 
coupling along the legs mixes ground sate and single¬ 
particle states with two-particle states and will give some 
corrections in the second order of perturbations, while 
SAE coupling on the rungs gives energy corrections al¬ 
ready in the first order on the coupling parameter: en¬ 
ergies of the ground state and of the S z = 0 compo¬ 
nent of triplet state are shifted (per dimer) by — jA rung , 
while energies of the S z = ±1 components are shifted by 
+\A rung . Zero field splitting of triplet levels appears, its 
type (easy axis or easy plane) depends on the coupling 
parameter A rung sign, which can be both positive and 
negative. 
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